SCALING LIMITS OF COUPLED CONTINUOUS TIME RANDOM 
WALKS AND RESIDUAL ORDER STATISTICS THROUGH 
- ( ^,-. MARKED POINT PROCESSES 

^ i ADAM BAECZYK AXD PETER KERN 

^ ■ Abstract. A continuous time random walk (CTRW) is a random walk in which 

, both spatial changes represented by jumps and waiting times between the jumps 

are random. The CTRW is coupled if a jump and its preceding or following waiting 
^v^j . time are dependent random variables (r.v.), respectively. The aim of this paper is 

to explain the occurrence of different limit processes for CTRWs with forward- or 
backward-coupling in Straka and Henry (2011) using marked point processes. We 
also establish a series representation for the different limits. The methods used also 
allow us to solve an open problem concerning residual order statistics by LePage 



<N 
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Oh 

(1981). 



1. Introduction 

Two i.i.d. sequences of M + -valued waiting times (J n )„ e N and of Revalued jumps 
(X„)„ e N yield two versions of a CTRW by 

od ' Nt Nt+1 

1 k=l k=l 



where N t := max{n G No : Ymc=i < t} is the number of jumps up to time t. The 
CTRW is coupled if the sequences ( J n )neN and (X n ) n£ N are dependent. Typically we 
assume, that the sequence ( J n , X n ) n6 N is i.i.d. with unknown dependence between the 
waiting time J n and the jump X n for fixed n G N. Using this dependence structure Sjv t 
is called backward-coupled CTRW whereas SV t +i is called forward-coupled CTRW. 
Both processes represent the position of a jumper at time t, but in the backward- 
coupled case the particle first waits for a time J\ before jumping to Xi, whereas in 
the forward-coupled case the particle jumps to Xi at time t — and then waits for a 
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time Ji and so on. CTRW processes were introduced in [22J to study random walks 
on a lattice and have been studied intensively over the past few decades. Today there 
is a wide field of possible applications for CTRWs. They are used in physics to model 
phenomena of anomalous diffusion [28J. The jumps can also represent movements of 
an ensemble of particles being transported over the earth surface in geophysics [27] 
or represent log-returns in finance [25J . A comprehensive study of limit theorems for 
coupled CTRWs has been initiated in [3] covering previously known special models 
[281 [TH [15] from physics. The limiting distributions of forward- and backward-coupled 
CTRWs have been investigated by Straka and Henry [29] using a continuous mapping 
approach on the space of their sample paths. A similar approach in a more general 
setting appears in [UJ. Straka and Henry prove that the limiting processes of coupled 
CTRWs in general differ when waiting times precede or follow jumps, respectively. 
Also the differences between the properties of these processes are not marginal, cf. 
[12] . Unfortunately neither the continuous mapping approach used in [29J nor the 
methods used in [TJ] are adequate to point out why different scaling limits occour. 
So a new approach to fill this gap is made using marked point processes here. Defining 
the time of the n-th jump by T n := ^22=1 ^ we ^ rs ^ study the limit behavior of the 
point processes which arise by marking each jump time with its occouring jump, 
respectively, i.e. we analyse 



It turns out that only the jumps with large norm contribute to the limit distributions 
of (11 .11) . as it is already known for real- valued partial sums which converge to an 
infinitely divisible r.v. without Gaussian part, cf. [1] and references therein. The 
methods used also solve an open problem concerning the convergence of residual 
order statistics by LePage, cf. [T7], [T8], [2S]- The scaling limits of the CTRWs 
can be determined by summing up the marks of the points in (11.11) which have a 
jump time occouring before time t. Hence in the scaling limit of forward-coupled 
CTRWs an additional big jump occurs compared to its backward-coupled version, 
which illuminates the difference between the processes. This approach also provides 
a series representation for the different scaling limits which might be of interest for 
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n 



(1.1) 




k=l 



k=l 
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simulation purposes. Since the resulting limit processes are not Levy processes, no 
efficient simulation algorithm is known yet. 



2. Preliminaries 

Let (J n ,X n ) ng N be an i.i.d. sequence of R + x Revalued r.v. Assume that (Ji,Xi) 
belongs to the generalized domain of attraction (GDOA) of a r.v. (D, A), where D 
is stable with index a G (0, 1) and A is full operator stable with index E G GL(M. d ) 
and without Gaussian part. Note that by Theorem 7.2.1 of [20] the real parts of the 
eigenvalues of E are greater than 1/2. By classic results [131 Theorem 14.14] and [21] 
Theorem 4.1] this implies that there exist regularly varying sequences (6„) ne N G RV\/ a 
and (A n ) ne ?q G RV-e such that the convergence 

/ \nt\ L«*J \ 

(2.1) b- 1 h, $^(^nX fc - E(A n X k l ]lAnXkl] < T )) A (D(t), A(t)) t > 



k=l k=l 



t>0 



holds in D([0, oo),M + x M. d ) for any r > such that for the Levy measure r] of A 
and the sphere = {x G M. d : \\x\\ = r} we have ?7(S^ _1 ) = 0, i.e. the sphere 

S^ _1 is a continuity set for rj. Here the process D(-) denotes an a-stable subordinator 
and A(-) denotes an operator Levy motion. Note that the drift term of the Levy 
process A depends on r. It is well known that we can choose r = oo if the real part 
of any eigenvalue of E belongs to (1/2, 1), since then E(Xi) exists. Moreover we can 
choose r = if the real part of any eigenvalue of E exceeds 1. Due to the spectral 
decomposition in [20], centering by truncated expectations in (12.1 ft is only necessary 
if some eigenvalue of the exponent E has real part equal to 1. 

We already stated, that only points with large norm contribute to the limit of the 
point processes in (11.11) . So we use a radial decomposition of the Levy measure 



1](A) = / / 1a(xv)t](cIx, v )da(v) 
igd- 1 Jo 

where a is a probabilty measure on the unit sphere of ~R d and (r/(-, v ))„ e §d-i is 
a weakly measurable family of Levy measures on (0, oo), cf. [24] . We also define the 
right-continuous inverse of rj{-,v) by 

(2.2) V*~i x i v ) := SU P{ M > : 7j([u, oo), v ) > x}. 
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With this notation we are able to give a series representation for the process A(-) in 
D([0, T] , R d ) for fixed T > by 



(2.3) lim V (ri*-(T l r fc , V fc )V fc l^ (r -i rfciVfc )> 6 ) -t / xdrj(x) 

£i ° \T-r k <t Je<\\x\\<r 



cf. [17] . [24] . where T n is the n-th partial sum of i.i.d. standard exponential r.v., 
( T n)neN denotes an i.i.d. sequence of uniformly U(0, l)-distributed r.v. and (V n ) ne N 
denotes an i.i.d. sequence with distribution a, with (r„)„ eN , (r n )„ gN and (V n ) neN 
being independent. 

Now it is well known that for a triangular array of infinitesimal row-wise indepen- 
dent M + -valued r.v. (Yk. n )i<k< n , nEN, converging to an infinitely divisible r.v. Y 
with associated Levy measure 0, only the extremes contribute to the limit distribu- 
tion, cf. [1] and references therein. This result coincides with the convergence 



(2.4) £%n,0 ^PRMW 



V . 
k=l 

in M p ([0, 1] x (0, oo]), where Pi?M(A\ ® 0) denotes a Poisson random measure with 
mean measure A\ <g> 0. Furthermore, it is well known that J2k&N £ (T k ,<t>^(r k )) * s a ^ so a 
representation of PRM{)\ <g> 0) in M p ([0, 1] x (0, oo]), where denotes the right- 
sided inverse of 0, cf. [21]. This fact can be understood by sorting the points on the 



left-hand side in @Q| 

n n 



n — fc + l:n 

fc=i fc=i v '" 



where (Yi m , . . . , F n:n ) denotes the order statistics of (Yi, n , . . . , with correspond- 
ing antirank vector (d\, . . . , d n ), i.e. the inverse permutation of the rank vector. 
Using Freedman's Lemma, cf. [10], one can easily verify, that the convergence 
(n _1 c4)fc e N — > ( r fe)fcGN holds in [0, 1] N , where d k = for k > n. Moreover, (2.4) 
in [Tj gives us F n _ fc+1:n — >■ 0^(r fc ). So the convergence in ( 12. 4p can also be es- 
tablished by analysing the convergence of the points (n _1 <i fc , F n _ fc+1:n ) 1 < fc < n . This 
approach can also be applied to the point processes in (11. ip . As the r.v. (X n ) nS N 
are Revalued one cannot use traditional order statistics. LePage suggested to use 
a normwise sorting, cf. [T7]. So for x\, . . . , x n G M. d we introduce the residual order 
statistics Xi :n , . . . ,x n]n by ||xi :n || < . . . < \[X n 



scaling limits of coupled continuous time random walks 
3. Convergence of residual order statistics 
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The convergence of the normalized residual order statistics A n X n _fe + i :n is still an 
open problem. LePage [T7J conjectures, that a generalisation of the one-dimensional 

case 



holds. As usual one sets Xfc :n = 0, whenever k < or k > n. This result has been 
proven in [18] for the case that the limit process A(-) is multivariate a-stable. In 
this case the right-sided inverse rj*~(x,v) defined in (12. 2ft is independent of v G S^ 1 
as the projection of the Levy measure rj is the same for every direction, cf. Theorem 
7.3.3 in |20j. Some years later a similar problem has been studied in [9] using a 
different norm || ■ \\h which respects the special structure of the operator E. In [26] 
the operator semistable case has been studied, but the result (13. ip also could only be 
established in the special case, that rf~(:r, v) is independent of v, which concides with 
the multivariate a-stable case. The author also supposed that the convergence (13.11) 
holds only in this case. We will show that the limit on the right-hand side in (13.11) 
has to be modified. The proof is based on the following lemma. 

Lemma 1. Let N n = X}fc=i £ x (n) ' n e ^ e a se Q uence of point processes in 

k 

M p ([—oo, oo] d \Kf) , the set of all point measures on the space [— oo, oo] d \K^ ; where 
:= {x G lR d : < e} denotes the compact e-ball in M. d . Suppose N n — Nq. If 
e < HX- "*!! < oo holds for every i G N almost surely (a.s.) and ||X^ || > ||X^ (w)|| 
for all 1 < i < j a.s. then the convergence 



Proof. The proof is based on a continuous mapping approach and Lemma 7.1 in [23] 




(3.1) 




v 




holds in (M d \K^) 




Define M C M p ([-oo, oo] d \K^) by 



p 



M : 



m : m = 




k=l 



Now we show that the mapping 



7T fc : M p ([-oo, oo] rf \Kf) -)■ [-oo,oo] d \K^ 



7TA: 
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£ xij ^ x P-k+i:P, xp-k+i-.p = for k < and k > P 



is continuous in m G M for every k G N. Let (m n ) n€N C M p ([—oo, oo] d \Kf) be a 
sequence of point measures converging vaguely to a point measure m = J2k=i e x {0) e 
M p ([— oo, oo] d \K^). Now choose n sufficiently large so that all points of m n lie inside 
of [—00, oo] d \IKg. By sorting the points of m n in descending order of their norm 

p 

II M\\ ^ ^11 (n)n . 

e ( n ), £ < II Xp II < ... < II || < 00 

i=l 

an application of Lemma 7.1 in [23] yields the convergence of the points 

(n) M\ , / (0) fO) 



^o.z,j IX] ,...,Xp 1 1 j-j ,...,Xp 

in (M d \K^) p . Now by the definition of the mapping Hk 

7Tfc(m n ) = 4. n) , 7r fe (m ) = 4. 0) for 1 < k < P 

TTfc(^n) = TTfc(^o) = for k > P 

holds and tt^ is continuous by (13.21) for every k G N. An easy application of the 
continuous mapping theorem yields the desired result. □ 

Lemma [1] allows to identify the distribution of the limit of properly normalized 
residual order statistics. 

Theorem 2. For k G N, T > and to G Cl define 
(3.3) d h {uj) = arg max ^(T^r^w), Vj(o;)) 

V^rfl^),...,^-!^) 

as £/ie argument of the k-th largest element of the set {^(T^r^u;), Vj(a;)), z G N}. 
TTien convergence of the residual order statistics 

(3-4) (A n X LnTJ _ fe+1:LnTJ ) fc6N A (r (T- 1 ^, V 4 )V 4 ) 

/jo/ds zn (M d ) N . 

Remark 3. Note that c4 is well-defined for all k G N, since the number of elements 
in the set {i G N : rf~(T'~ :I Tj, Vj) > e} is finite a.s. for all e > 0. Moreover 
theorem [2] does not contradict any of the results proven in [9], [18], [26 J. If A(-) is a 
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multivariate a-stable Levy process, the monotonicity of the mapping x h- > tj*~(x,v) 
yields dk{w) = k a.s. for all k G N. 

Proof of Theorem First we have to determine the limit of the truncated point pro- 
cess 

\nT\ 

£(A„X fc l||_4 nXfc ||> e ); 

fc=l 

where e > has to be choosen such that ^(S^ -1 ) = holds. By Theorem 3.2.2 in [20] 
assumption (12.11) yields the vage convergence 

(3.5) [nT\F(A n X k e-) ^T- V (-) 
in M d \{0}. Hence the convergence of the point processes 

\nT\ 

fc=i 

holds in M p ([-oo, cxf \{0}). Now by [24] J2keN^(T^r k ,v k )v k = PRM(T ■ r,). Ap- 
plying the a.s. continuous restriction functional 

(3.6) tt' : M p ([-oo, oo] d \{0}) -> M p ([-oo, oo] d \K^), m ^ m m)C , 
the continuous mapping theorem yields 

\nT\ 



( 3 - 7 ) £ ^» x fc 1 iiA„x fc n> e y^g^(T-ir fc ,v fc )v. l. 



fc=i fceN 
The continuity of 7r' is proven in [5] for instance. Now the points of the point process 
on the right-hand side of (13.71) have to be ordered in descending order of their norm 



fceN fceN 
An application of Lemma Q] yields the convergence of the points 



A n X y nT j _ k+1 . LnT j 1 ||A„X LnTJ _ k+1: LnTJ || >e 



fceN 



V 

-(T—I-Ft |Vj )>£ 



fceN 

in (M a! \]K^) N . The desired result follows by taking the limit as e 4- and an easy 
application of Theorem 4.2 in [1]. □ 
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4. Convergence of associated point processes 



Now that the limit distribution of normalized residual order statistics is identified 
we can study the associated marked point processes 

LnTJ [nT\ 



k=l k=l 

In the uncoupled case convergence results for this processes can be established with 
a continuous mapping approach using the time deformation defined in [231 (8.29)]. 
Since the continuity of this time deformation demands the processes A(-) and D(-) 
to have a.s. no common jumps, which is not necessarily fulfilled in the coupled case, 
this standard methods connot be applied in our case. So we use a sorting argument 
like in (1231). 



Lemma 4. Let (r n ) neN; (T n ) ngN and (V n ) neN be as in ( \2.3\i . Then the convergence 
of the associated point processes 

e {b~ 1 T k ,A^ k ) y ^2 £ (D(T-r k ),^(T^r k ,V k )V k ) 

k=l keN 
LnTJ 

( 4 - 2 ) £ {K 1 T h - lt A n X k ) ► £ (D(T-T k -),^(T-iT k ,V k )V k ) 



k=l keN 



holds in M p ([0, oo) x [— oo, oo] d \{0}) for every T > 0, where D(x-) denotes the 
left-hand limit of the process D(-) in x. 

Proof. Choose T > arbitrary. We start by sorting the points of the associated point 
process 

\nT\ \ (\nT\ 

Y e (b- 1 T k ,A^ k ) = I E £ (VT d LnTJ _ fc+1 ,A„X LnTJ _ fc+1:LnTJ N 
k=l / \ fc=l 

Again (di, . . . , d^ nT ^) denotes the antirank vector of the r.v. (Xi, . . . , X|_ n yj). The 
normalized residual order statistics A n X„_ fc+ i : i nT j have already been studied in The- 
orem [21 It remains to determine the limit distribution of 



d-lnT] -fc+1 



\nT\ ■ "L nT J -fc+i 



1=1 1=1 
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Since X 1; . . . , X|_ n rj are i.i.d., (n _1 <i fc ) fceN — > {T-Tk)kes suggests that the convergence 
( 4 - 3 ) (K lT d lnTi _ k+1 )ken {D{T ■ r k )) km 

holds in [0, oo) N . But this result cannot be established with a traditional continuous 
mapping approach. The mapping ir t : D([0, oo),R d ) x IR + — > R, (x,t) i— > x(t) is only 
a.s. continuous if x is a.s. continuous in t. Also classical transfer theorems, cf. [7], 
[8], are not helpful because they require independence of the summands and their 
quantity or a stochastic convergence of the normalized antirank vector, cf. [2]. Since 
none of these conditions is fulfilled another approach is used. 

Let (di, . . . , d\ n T\) denote the associated antirank vector of the waiting times 
(Ji, . . . , J\nT\)- Since the joint convergence of the well-centered and normalized se- 
quential partial sums to the process (£)(•), A(-)) holds, one can easily prove conver- 
gence of the normalized antirank vector 

(n 1 (d[ nT j_ k+ i) ke ?q,n 1 {di n T\-k+i)km) — * ((T ' r i)ieN) (T ■ Tj)^) 

where (r n ) ne N, (T n ) n m denote two i.i.d. sequence of U(0, redistributed r.v.. Note 
that the sequences (r n ) n€ ^ and (r„)„ gN are not independent in the coupled case. As 
a consequence of the convergence of the antirank vector, the convergence of indicator 
functions 

(l ~ ) V ) (1 ) 

\ n 1 d lnT} _ j+1 <n 1 d lnT1 _ i+1 J V T i< T i/j e ® 

holds. An application of Basu's lemma, cf. [161 Theorem 5.1.2], yields the indepen- 
dence 

(4.4) Mdi, . . . d\nT\) , (di, . . . d[ n T\)} -L (Jl : [nTJ, • • • , J[nT\:[nT\) 

for every fixed n G N, which proves that the joint convergence 



[nT\-j+l:[nT\i-d lnn _ J+1 <d VnT ^ l+1 ) jm 



holds, where {f] a )~ l is the right-sided inverse of the Levy measure associated with 
-D(l) and (r n ) n6N denotes a distributional copy of the sequence (r n )„ eN . Note that 
the sequences (r„) neN and (r„) neN are also not independent in the coupled case. Now 
summation verifies 

(4-5) b n lj lnT\ -l+l: [nTj <n^d LnTi 

zeN 
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\nT\ 



E b n lj \ 



1 = 1 



LnTj 



."'J-'-l:L"TJ J -n-id LnT j_, +1 <n-id LBT j_ m 

V 



(4.6) =b- i j2 J i l i<^ Ti ^ = K 1 £ Ji^EOT" 1 ^') 1 ^^- 

z=i i=i zeN 

Using the Ferguson-Klass series representation of the process D(-), cf. [6], one iden- 
tifies the right-hand side in ( 14. 6 p as series representation of D(T ■ Tj). Since the 
convergence of antiranks holds simultaneously, we have proven 

/ [n . VTM +l j \ 



(4.7) 



K 1 £ J < 



V 



1=1 



V 



(D(T ■ Ti)) i 



/ 



■ieN 



111 



} . Again by the joint convergence of the properly normalized and scaled se- 



quential partial sums to the process (£?(•), A(-)), the independence (14.41) and Theorem 
121 convergence of the points 

holds. Since the mapping x i— > e x is continuous, the convergence of points yields the 
convergence of point processes 

v 



' /b " lTd LnTj-i + l' A " X L"rj- l+ l:LnTj y y _ 



D(r^,),^(r- 1 r ? .,v 5 -.)v ? 



in M p ([0, oo) x [— oo, oo] a! \{0}). Now summation and an easy application of theorem 
4.2 in [4] yields 



LnTj 

E 



£ ( bn lTd lnTi -k+1 ' A « X L"TJ -k+1: LnTj y 



(D(T.r fc ),^(r-i r(rfc ,V 4 )V 4 ) 



fe = l L J I , 

in M p ([0, oo) x [— oo, oo] d \{0}). Now we need to reverse the order of the points again. 
We introduce 



as inverse of dy.. Note that is well-defined for all k G N by Remark [31 Moreover, 
since (r n )„ eN and (r n )„ eN are independent, an easy application of the desintegration 
formula shows, that (j rn ) n eN is also i.i.d. and U(0, 1) distributed. Since the conver- 
gence in ( 12. 41) towards (12. 3ft can be proven with the same technique, the sequences 



SCALING LIMITS OF COUPLED CONTINUOUS TIME RANDOM WALKS 



11 



( r rn)neN, (r n ) n£ N and (V„) ne pj are also independent. So we define r n := r rn for all 
n G N and obtain 



[nT\ 

E 

k=l 



£ (b- 1 T k ,AnX k ) 



\nT\ 

E 



k=l 



" ' K lTd [nTi -fc + 1: LnTj ' A " X L«Tj -fc + 1: |nT| 



^ £ (.D[T. Vt ),rj<-(T- 1 r k ,V k )V k )- 

ken 



Hence we have proven (14.11) . In order to prove fl4.2[) the limit distribution of 



h~ lr T 

u n - t d[„TJ-i+l- 1 



" Ln rj-i+i_ 1 

n 

k 1 E Jfc = 6 « 1 E 



fe=i fc=i 

has to be analized. Denoting G([0, oo),M, d ) the space of all left-continuous functions 
with right-hand limits from [0, oo) to M. d one easily proves that the mapping T : 
D[0, oo) — > G[0,oo), I—)- x(t— ) is Lipschitz-continuous with Lipschitz-constant 
one and hence continuous. Since 

/ LntJ \ fn*l-l 

r E J * = ^ Jfc 

yfc=i y fc=i 
holds, the continuity of T suggests, that the convergence 

(4-8) | b-' E 



k=l 



V 



(D(T-n-)) 



J, 



holds. But since (" 14. T[) could not be proven with the continuous mapping theorem, we 
use the above arguments again to obtain (14. 8p . We note that 



-1 



LnTj 



V E Jl = J2 b n lj \ 



1=1 



1=1 



[nT\ 



L" T J-'+l:L"*J J -d LnTJ _j +1 <d L „ TJ _ i+1 -l 



(4.9) = E ^Ann-i+i^K-^^Kn-^^ — > E^'( T ^Ir-nKT-n 
i=i len 

holds. Again we identify the limit on the right-hand side in (14. 9p as a series repre- 
sentation of D(T ■ Tj— ). As already stated this yields the desired result (14. 2p . which 
completes the proof. □ 
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5. Scaling limits of coupled CTRWs 

In this section we are now able to identify the scaling limits of coupled CTRWs 
using the limit theorems for their associated point processes stated in Lemma HI We 
need to introduce the set 

S :={T ER + : ¥(D(T ■ r<) = T) = for all i E N} 

for technical reasons. Due to selfsimilarity, one can easily show that the equality 

D(xt ■ n) = x 1/a D{t ■ n) 

holds for all % e N and x, t £ R + . Hence the set S is dense in R + . 



Theorem 5. Let E(t) := inf{a; > : D(x) > t} denote the hitting-time process 
associated with D(-). Then convergence of the backward- and forward- coupled CTRW 

^(A n X fc - E(A n X fe l||A„X fe ||<r)) 
k=l 

(5.1) 



v 



lim £ (fr^.V^V.VtT-r.v^s) - E(t) [ 

El \D(T-r k )<t 



x drj(x) 



and 

J2 (AA - E(A n X fc l||A n X fe ||<r)) 

fe=i 
(5.2) 



v 



\D(Tr k -)<t Je ^ x 



x dr](x) 



holds in D([0, T], R ) for every T £ S and every r > such that r/CBf ) = 0. 

Proof. Choose T e S arbitrary. Similar to (13. 6p we define another a.s. continuous 
restriction functional 

n : M p ([0, oo) x [-oo, oo] d \{0}) ->• M p ([0, oo) x [-oo, oo] d \Kf), 

7T (m) := m,|[0, oo )x[-oo,oo]' 1 \Kf 
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for e > such that ^(S^ 1 ) = 0. Moreover we define the summation functional 
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X : M p ([0,oo) x [-oo,oo] d \K^) -> D([0,T] 

x X! e M ) (*) = 

VfceN / \t fc <t / 



te[o,T] 



which is a.s. continuous in the point 



KD(T.r fc ),f?-(T-lr fc ,V fc )V fc l^ (T -i rfc , Vfc) >J 

for every T G S. A proof of the continuity of \ is given in [231 Sec. 7.2.3] for the case 
d = 1 and can easily be modified to hold for d > 1. So we apply the a.s. continuous 
mapping x 07r to the associated point processes in Lemma HI Considering the equality 
{T n <t} = {N t > n} we receive 

* 07r \J2 £ (b^T k ,AnX k ) ) (*)= ( E ^ X ^HA„X fc ||> £ 



te[o,T] 



A n x k i llAM > e \ — >x°k{J2 £ 

- fc=1 Jte[o,T) ^ keN 

(5.3) = ^ ^(T _1 r fc , V fc )V fc l^( T -i rjk ,Vfc)> e - 

D(T-Tk)<t 

To study the centering constants we use the convergence 



(D(T-r fe ),^(T-ir fe ,V fc )V ft ) 



n 7 t>0 



t>0 



in D([0, oo),R + ), proved in Corollary 3.4 of [21]. Considering (13. 5 j) this yields 



fc=i 



VE(A„Xa [£ , T ](P n X fc ||)) = / x ndF 4 "* 1 ^)) 

£7(f) / x rfr/(x) J 



'e<|[a;]|<r / t>0 
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in D([0, oo), IR + ). Since the process E(-) has a.s. continuous sample paths, Theorem 
4.1 in [30] allows us to put this and ( 15. 3 1) together. We obtain 



^ (^nX fc l||A n X fc ||>e - E (v4 n X A .l e <|| AnXfc ||<r)) 



k=l 



te[o,T] 



V 



5^ {v^(T k , V fc )V fc l^ < - (T -ir feiVfc )> e ) - E(t) / x dr](x) 

\D(T-r k )<t Je<\\x\\<r 



te[o,T] 



in -D([0, T], P d ). Taking limits as e \. this yields the desired result ( 15. ip . By Theorem 
4.2 of [4j it remains to show 



limlimsup PI sup 

e i° n-Hx> \ 0<t<T 

N tbn 



2_j { A nKk1\\A n y. k \\>s - E {A n X. k l £ <\\ An x k \\< T )) 



k=l 



~ OnX fc - E (^ n X fe l||A„X fc ||<r) •) 



k=l 



> 5 = 



for all 5 > 0. Using a version of the Kolmogorov-inequality for integrable stopping 
times given in the Appendix and the norm-inequality || • || < || • ||i we obtain 



P( sup 

0<t<T 
N tb„ 



N tbrl 



22 (^nX fc l|| AnXfc ||> £ - E (AiX fc l £ <|| AnXfc ||< r )) 



k=l 



(4»Xfc - E (A n X fc l||A n X fc ||<r)) 



k=l 



> 5 



( 



< P max 

\l<j<N T b n 



> 5 



< P V max 

» -f-* l<j<N Tb 
,i=l 

d 

< P [ I J max 



y j A n X fc l|| AwXfc || <£ - E (AiXfeln^Xfel^e) 

^(A n X fc )«l|| A , iXfc || <£ - E ((A n X fc )«l|| AnXfc || <£ ) 
fc=i 

^(A n X fc )«l|| AnXfc || <£ - E ((A n X fc )«l|| AnXfc || <£ ) 



fc=i 



> s 

6 



> 



d 

i=l 



max 

l<j<N Tbn +l 



J2(A n x k fh UM<£ - E {(A n x k fh UM<e ) 



k=l 



d 



> 



d 
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-2 d 



< -. 



<'5 



E(N Tbn + l)j^Var {{A n X x fh UM<£ ) 

- ' ' i=l 

! d 

E(iV T6n + l^E (((A n X 1 Y% An ^ ll]<E ) 2 



i=i 

where x^ denotes the n-th coordinate of the vector x. Note that we have to take 
N t b n + 1 since N t b n does not fulfill the conditions of Lemma [71 Now Theorem 9 in [19] 
states that E (N t + 1) can asymptotically be expressed by the integrated tail of the 
distribution function of J\ 

E(iV t + l)r(2-a)r(l + a) / 



£{1 - Fj&Vdx' 

where T denotes the gamma-function. By Karamata's Theorem the limit behavior of 
the function 1 1— > J *(l — Fj 1 (s))ds can be expressed by Fj 1 

Jo 1 - a 

Putting this together we obtain 

E (N t + 1) ~ (1 - Fj.it))- 1 • r(l - a)' 1 • r(l + a)" 1 . 

Defining C := (r(l — a) T(l + a) r] a ((T, oo)))^ 1 for abbreviation, the inequality 
\x^\ < yields 

-2 



- j limlimsupE (N Tbn + 1) VE ( {{A n ^H UM<e \ 
' i=i 
/fi\~ 2 d 

< limlimsup I - C- limlimsup nE ( ( (A^Xi ) ^> ) 2 1 , (AtiXi } (i) , <£ ) . 

e +" n— >oo \" J ■ , e ^ n— >oo 

x ' i=l 



So it remains to show 

limlimsupnE(((A n X 1 )«) 2 l |(AnXl)(l )| <£ ) = 

n— >-oo 

for all 1 < i < d. Defining the tail and truncated second moment of Xi in direction 

v by 

V(r,v) :=F(\(X 1 ,v)\>r), U(r,v) := E((X 1 ^) 2 l KXl ,,)|< r ) 
for every v G and every r > an easy calculation shows that 
E (((A n X0 w ) 2 l| ( ^ Xl) «|< £ ) = r 2 n U(r^e,v n ) 
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holds. Here r n > and v n G S d_1 are taken such that A* n ei = r n v n holds for every 
n G N, where A* is the adjoint of A and ei, . . . , e<f denotes the standard basis of lR d . 
With this notation we have to analyse 

5.4 n-r n E7 [r n e,v n ) = e — — — - — • — — — — • n ■ V[r n ,v n ). 

One easily verifies that the third factor in f!5 .4jl is bounded by r]({x G W 1 : \x^'\ > 1}). 
Since the real parts a\ < . . . < aa of all eigenvalues of the operator E are greater 
than 1/2 we can find e > such that 2 — e — a^ 1 > holds. So an application of 
Theorem 6.3.4 of [20] yields the existence of a constant C\, such that for neN large 
enough the second factor in (I5.4|) is bounded by C\E 1 ~ e ~ a ^ . Finally, Corollary 6.3.9 
in [2U] yields that the first factor in (15.41) is bounded by a constant C%. Putting things 
together this completes the proof. 

The proof of the convergence (15. 2p works the same way. □ 

Theorem [5] provides a series representation for the limit distribution. This repre- 
sentation might be useful for simulation purposes. Now it is of interest to identify the 
scaling limits with the ones stated in [29]. The arguments are based on the following 
equalities: 

(5.5) {D(x) <t} = {x< E{t)}, {D(x-) <t} = {x< E(t)}. 

The left-hand side of (15.51) is already proven in (3.2) of [21J. For the proof of the 
right-hand side assume D(x—) < t holds. So we have D(y) < t for all y < x. Hence 
x < E(t). Otherwise if D{x—) > t holds, there exists an e > such that D(y) > t 
holds for all y > x — e. Hence E(t) < 

Corollary 6. Let A(t— ) + denote the right- continuous version of the process A(t—), 
i.e. A(t—) + is an element of D[0,oo). Then the convergence 

N tb n 

(5.6) ^(A„X, - E(^ n Xa||A„x fc ||<r)) ^ ME(t)-) + 

k=l 

and 

(5.7) Yl ( A " X x ~ ^A.Xfelii^n^)) A A(E(t)) 

k=l 

holds in D([0, oo),M d ) for any r > such that ^(S^ 1 ) = 0. 
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Proof. The convergence (15.71) can easily be verified applying the right-hand side of 
(15.51) to the series representation ( 12.31) . To prove (15.60 we apply the left-hand side of 
(1531) to flSHJ) and obtain 

Y,(A n X k -E T (A n X k )) 



k=l 
V 



lim {^(T-'T,, V k )V k l^ { T-ir k ,v k )>e) ~ E(t) [ 

£i ° \D(T-T k )<t Je<\\x\\<r 

= E {v^iT^Tk, V fc )V fc l^ (T -ir fclVfc )>J 

\T-r k <E(t) 

(5.8) 

+ ^ (^"(r- 1 r fc ,v*)v fe V(T-ir fc ,vo> E ) - E V) [ xd v(x) 

t - \—t •/e<||a;||<T 



x dr](x) 



D{T-r k )=t 

As we already stated 



^ \T^<B(t) ^<||*||<t 



x dr](x) 



is a series representation of A(E(t)—). The extra summands only have to be consid- 
ered if a jump occurs at a time t with D(T ■ r k ) = t. This yields the right-continuity 
of the limit and we have proven (15.61) . □ 

Appendix A. A generalization of Kolmogorov's inequality 

The following generalisation of Kolmogorov's inequality can be shown by standard 
techniques. However, we were not able to find a suitable proof in the literature and 
will only give a sketch of proof. 

Lemma 7. Let (Y n ) n ^ be i.i.d. with E(Yi) = and T be an No-valued integrable 
stopping time with respect to the filtration JF n := o~(Yi, . . . , Y n ). Then 

k 

>8) < <T 2 • E(T) • Var(Y y 



P I max 

Kk<T 



En 



holds for alle 5 > 0. 
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Proof. First we restrict our attention to the truncated stopping time TAn. An easy 
calculation shows that 



(TAn \ 
g (S k - 5 fc _x) 2 J = E(fi£ AB ) 



holds, where £„ := ^fc=i denotes the n-th partial sum. Defining 
M := 0, M k+1 :-. 



S k +i, if maxi<,< fc \ Sj\ < 5 
Mk, else 



the same calculation shows 

E (M k - il4_0 2 ) = E(M| A J - 2 ■ E ( J^iMk ~ M^S^lkKTAn 
\fc=l / VfceN 

Using the definition of M n one can show 



^TAn) 



Now |Mfe — Mfc_i| < | jSfc — holds for all k G N. Hence the Markov-inequality 

yields 



P 



( max 1^1 > ^ = P (| -Mr An I > 5) < <T 2 E (M 2 An ) 

/TAn \ /TAn 

-8- 2 E[J2 (Mk - M fc _ x ) 2 <r 2 £g (5Jk - S fc -i) 2 ) = <T 2 E (5« Afl ) 



vfc = l / \fc=l 

An application of Wald's inequality yields the desired result for TAn. The generali- 
sation for the stopping time T follows by the martingale convergence theorem. □ 

Remark 8. Let (Z n ) ne N be a sequence such that (Y n , Z n ) n ^ is i.i.d.. Then Lemma [7] 
also holds replacing T n by := o-((Y u Zi), . . . , (Y„, Z n )). 
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